The characteristics of neutrinonuclear reactions 

at = 1 - 3 GeV 



SKAT Collaboration 

N.M. Agababyan-'^ , V.V. Ammosov^, M. Atayan^, 
N. Grigoryan^, H. Gulkanyan^, A. A. Ivanilov^, 
Zh. Karamyan^, V.A. Korotkov^ 



^ Joint Institute for Nuclear Research, Dubna, Russia 
^ Institute for High Energy Physics, Protvino, Russia 
^ Yerevan Physics Institute, Armenia 



YEREVAN 2004 



Abstract 



For the first time, the characteristics of the charged current neutrinonuclear interactions 
are investigated at = 1 - 3 GeV using the data obtained with SKAT propane-freon 
bubble chamber irradiated in the neutrino beam at Serpukhov accelerator. The E^, 
- dependence of the mean multiplicities of different types of secondary particles and 
their multiplicity, momentum and angular distributions are measured. 



1 Introduction 



The experimental data on the characteristics of neutrinonuclear reactions at <3 GeV 
are of interest, in particular, in view of current and anticipated experiments on the neutrino 
oscillations (see, e.g., P). The detailed data are available only in the low energy range of 
Ejj <1 GeV, where the quasielastic reaction vn fi~p on the nuclear neutrons dominates 
The aim of this work is to obtain, for the first time, detailed experimental data on 
neutrinonuclear reactions in the range of 1< E,j <3 GeV. Section 2 is devoted to the exper- 
imental procedure of the selection of the neutrino interaction events and the reconstruction 
of the neutrino energy. The experimental results are presented in Section 3 and summarized 
in Section 4. 



2 Experimental procedure 

a) The event selection criteria 

The experiment was performed with SKAT bubble chamber |3] , exposed to a wideband neu- 
trino beam obtained with a 70 GeV primary protons from the Serpukhov accelerator. The 
chamber was filled with a propan-freon mixture containing 87 vol% propane (C^Hs) and 
13 vol% freon (CFsBr) with the percentage of nuclei H:C:F:Br = 67.9:26.8:4.0:1.3 %. The 
nuclear interaction and radiation lengths of the mixture were A/ = 149 cm and Xq = 50 
cm. The volume of the chamber was 6.5 m^, the chosen fiducal volume 1.73 m^. A 20 kG 
uniform magnetic field was provided within the operating chamber volume. 
Charged current interactions, containing a negative muon were selected. A muon was identi- 
fied as a particle that possessed the highest transverse momentum among negative particles 
and did not suffer a secondary interaction in the chamber. Other negatively charged parti- 
cles were considered to be vr" mesons. The overwhelming part of protons with momentum 
below 0.6 GeV/c and a fraction of those with momentum up to 0.85 GeV/c were identified 
by their stopping in the chamber. Other positively charged particles were assumed to be 7r+ 
mesons, except for cases explained below. 

In order to reduce the uncertainty in reconstructing the neutrino energy E^, the events were 
selected for which errors in measuring the momenta of all charged secondaries and converted 
photons were less than 27% and 100%, respectively. The mean relative error < Ap/p > in 
the momentum measurement for muons, protons, pions and gammas was, respectively, 4%, 
6%, 10% and 21%. Each selected event was assigned a weight that took into account loss of 
events. The mean weight of the sample used in the present study was 1.18. 
To diminish the background induced by neutral and charged particles, the following require- 
ments were imposed on selected events: 

- The net charge of secondary hadrons is positive. 

- The visual energy of neutrino E"™** = E^j^ + u^s > 0.6 GeV, where i?^ is the muon energy, 
Uyis is the visual energy transferred to hadrons (including decay photons and neutral strange 
particles, if any). 

- The value of the summary momentum of detected secondary particles was | \> 0.7 
GeV/c, while its longitudinal component pi (with respect to the neutrino direction) exceeded 
the transverse component pr (the latter coinciding with the summary transverse momentum 
of undetected secondary particles). 
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- The minimal neutrino energy defined from tlie kinematics of the quasielastic reaction 
un — >• fj,~p, Ef = {mjqE^ — Q.bm^) / {mjq — + p^cos-d^), did not exceed more than by 
30% [2] the upper hmit (3 GeV) of the considered range of E^: Ef < 3.9 GeV. Here m^v 
and rUfj, are the proton and muon masses, and -i?^ are the muon momentum and angle 
with respect to the neutrino. 

The distribution on i?"** for selected events in the range 0.6 < E^^^ < 4 GeV is shown in Fig. 
la. The measured value of i?™* is, in general, the minimal value of the true neutrino energy, 
and hence to be corrected for the contribution from undetected gammas and neutrons. An 
exception are the events satisfying the kinematics of exclusive reactions without secondary 
neutral particles (see the next subsection). For a part of these events, E^, can be even smaller 
than provided that the proton hypothesis for an unidentified positive particle matches 
better the kinematics of an exclusive channel, than the 7r+ hypothesis. 



b ) The identification of exclusive channel 

In order to avoid a unnecessary energy correction for undetected secondaries, each event was 
analysed for its belonging to the following exclusive channels without neutral secondaries: 



un ^ p 

up fl^p TT^ 



(1) 
(2) 
(3) 



with an arbitrary number of accompanying identified protons (originating from the nuclear 
cascading). In (l)-(3), p (vr^) stands for an identified proton (vr^ meson) or a positively 
charged particle to be identified as a proton or -n^ meson from the kinematical fit. 
Figs. 2a and 2b show the distributions on the azimutal angle between vectors "p*^ and 
= in the transverse plane and on the squared transverse momentum p"^ for 

events topologically consistent with reactions (l)-(3). The observed peaks near ~ 180° 
and ~ indicate on a significant fraction of events satisfying the kinematics of reactions 
(l)-(3). For further analysis the events with (p > 155° andp"^ < 0.15 (GeV/c)^ were selected, 
the chosen boundaries roughly corresponding to violation of the kinematics of reactions (1)- 
(3) due to the Fermi-motion of the target nucleon and the measurement errors. For these 
events, the distribution on the dimensionless variable = [E^^'^ — p^) / E'^^'^ , characterizing 
the longitudinal momentum-energy disbalance of the reaction, is plotted in Fig. 2c. As it 
is seen, the distribution is strongly peaked at ^ 0. The events with | |< 0.15 were 
assumed to be belonging to the exclusive reactions (l)-(3). The number of these events is 
126 (the weighted number 140.2), the distribution of which on E^, is plotted in Fig. 2d and, 
separately for channels (1), (2) and (3), in Fig. 2e. The curves in this figure are obtained 
by approximation of available experimental data for up and un interactions j3] (taking into 
account the proton and neutron contents of the target nuclei) convoluted with the expected 
neutrino flux spectrum on the bubble chamber SKAT [S]. The normalization of these curves 
is chosen to fit our data on the main exclusive channel (1). As a result, the predictions 
for channels (2) and (3) overestimate the data, probably due to the intranuclear interaction 
effects for pions produced in the primary uN interaction (see Section 3 below). 
The distribution on the invariant mass W of the hadronic system for exclusive channels (1) 
and (2) is shown in Fig. 2f. The clear peak at W ^ ttin, with a typical width ~ 0.15 GeV 
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(caused by Fermi- motion and measurement errors), corresponds to the quasielastic reaction 
(1). The W- distribution for the reaction (2) reflects probably the production of the 
states. 

c) The reconstruction of the neutrino spectrum 

Only a small part (about 28%) of the total event sample can be identified as exclusive 
channels (l)-(3). The rest part contains undetected secondaries, for which a correction has 
to be introduced to estimate E^. The mean value u of the corrected energy transferred 
to hadrons (at given narrow intervals of Uyis) is approximated as z/ = a + buyis, where the 
coefficients a and b are found by means of the procedure applied in [6j a = 0.15 ±0.24 GeV, 
b = 1.07 ± 0.05 at E^J' > 2 GeV and a = 0.5 ± 0.3 GeV, 6 = 1.5 ± 0.4 at E^J' < 2 GeV. 
The reconstructed spectrum of Ei, in the range 1-3 GeV is presented in Fig. lb, while 
Fig. Ic shows the spectrum of the beam neutrino evaluated using the total cross sections 
[7] of reactions up —>■ fi~X and un —>■ fi~X and taking into account the proton and neutron 
content of the target nuclei. As it is seen, the latter spectrum is consistent with the expected 
neutrino flux spectrum (at E^ > 1.5 GeV) on the bubble chamber SKAT [S]. 

3 Experimental results 

a) The E^- dependence of the mean characteristics of produced particles 

Fig. 3 presents the E,y - dependence of mean multiplicities of secondary particles: protons, 
TT"^ mesons, 7r~ mesons, charged pions, 7- quanta, mesons and neutral strange particles 
(V^). The mean value of < n^o > is estimated as < Uj^o >= 0.5 < > w^, where is the 
mean weight of detected 7's deflned by their conversion probability in the chamber. As it is 
seen, all multiplicities increase with E^, except the mean multiplicity of identifled protons, 
probably due to the fact that, with increasing E^, a larger fraction of recoil protons acquires 
momenta Pp > 0.6 GeV/c and hence can be detected as a unidentifled positive particle. The 
contribution of these protons is partly taken into account in Fig. 3a, where the full circles 
stand for < rip > including also 'kinematically' identified protons from reactions (l)-(3). The 
latter will be included below in the all data concerning protons. The solid curves in Figs. 
3b, 3c and 3f are the approximation of the experimental data jl] on the mean multiplicities 
of 7r+, vr" and vr" mesons in uN interactions. It is seen, that the values of < n7r+ > at 
E^ = 2.3 - 3 GeV exceed those in uN interactions, probably due to a significant (~ 30%) 
contamination from the unidentified protons. On the other hand, the values of < ?t,^+ > 
at i?,^ = 1 - 1.9 GeV and < ra^- > and < n^o > in the whole Ei, range are reduced as 
compared to those in uN interactions. This reduction can be, at least partly, attributed to 
the absorption of pions on a pair of nucleons within the nucleus, nlNN) —>■ [NN'). The 
dashed curves in Figs. 3b, 3c and 3f are the results of calculations incorporating this process 
(see for details j2] and references therein). As it is seen, the predictions are in reasonably 
agreement with the data on < > at -E;^ = 1 - 1.9 GeV, as well as on < Ut^- > and 
< n^o > in the whole energy interval, especially if one takes into account the uncertainties 
(reaching ± 20%) in the approximation of the experimental data on the mean multiplicity 
of pions in uN interactions. 

The E^ - dependence of the mean momentum of different types of particles is presented 
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in Fig. 4. This dependence is more expressed for charged pions and especially for muons, 
being rather moderate for protons and gammas. Fig. 5 shows the E^, - dependence of the 
mean fraction /j of the neutrino energy transferred to the muon (/^ = E^/Ey)^ protons 
{fp — YjTp/Ey), charged pions (/^ = J2Et,/Ei,), detected gammas (/^ = Y^E^/E^,) and 
undetected neutral particles {fmiss — {j^ — J^vis) / E^) ■ It is seen, that depends on E^, rather 
weekly, being enclosed between 0.41 < < 0.47. The proton contribution smoothly de- 
creases from fp = 0.13 to 0.05, while that for charged pions increases from f^^ = 0.12 to 0.32. 
The contribution of the detected gammas is very small, amounting = 0.01 - 0.04, and 
increases with E^,. The energy fraction carried by the undetected neutral particles decreases 
from fmiss ~ 0.3 up to ~ 0.1 with increasing E^. 

The mean values of the invariant mass W of the hadronic system and the squared transfer 
momentum are plotted in Fig. 6. A significant rise of < > is observed, while the E^,- 
dependence of < > is rather weak. 

b) The multiplicity, momentum and angular distributions of produced particles 

Figs. 7a - 7f present the multiplicity and momentum distributions for protons, as well as 
separately for those produced in the forward {cos6p > 0) and backward {cosOp <0) direc- 
tions with respect to the neutrino beam. These distributions, as well as the proton angular 
distribution (Fig. 7g), are presented for the whole range of E',^ = 1 - 3 GeV, because the 
characteristics of protons were found to not depend noticeably on E^,. The multiplicity dis- 
tributions, shown in Figs. 7a - 7c, arc compared with the poissonian distributions calculated 
at the measured values of < rip >= 1.17 ± 0.05, < rip >= 0.83 ± 0.04 (at cos6p > 0) and 
< rip >= 0.33 ± 0.02 (at cos9p < 0). As it is seen, the tail of experimental distributions de- 
viates from the poissonian one. The momentum distributions (Figs. 7d - 7f) can be roughly 
approximated by the dependence ~ exp{—ap^) with a = 3.5 ± 0.2 (at all cosOp), 2.8 ± 0.2 
(at cosOp > 0) and 4.6 ± 0.4 (GeV/c)"^ (at cos9p < 0). The distribution on cos9p (Fig. 7g) 
exhibits a linear dependence ~ 1 + j3cos6p with a slope parameter f3 = 0.82 ± 0.05. 
Figs. 8 and 9 show the momentum and angular distributions of muons and charged pions 
in two intervals oi E^ — 1 - 2.5 and 2.5-3 GeV with approximately equal statistics. The 
momentum distribution of muons (Figs. 8a and 8b) is peaked at ~ 0.5 GeV/c and falls 
steeply with increasing p^ for the range E^ = 1 - 2.5 GeV, being, however, much flatter for 
the range Ei^ = 2.5 - 3 GcV. The angular distribution of muons is similar for both Ei^ ranges, 
being strongly peaked at cos6'^ ~ 1. The momentum distribution of pions is also peaked at 
Ptt ~ 0.5 GeV/c and falls (more steeply for the range E^ = 1 - 2.5 GeV) with increasing p^^. 
Their angular distribution is peaked at cos^^ ~ 1 (more strongly for the range E',^ = 2.5 - 3 
GeV). 

The momentum and angular distributions of 7- quanta, both corrected and not corrected for 
their detection efficiency, are shown in Fig. 10, but now, due to the lack of their statistics, 
for the whole range of E',^ = 1 - 3 GeV. The momentum distribution has a maximum at 
Pj ~ 0.15 GeV/c (more than twice smaller than that for charged pions), while the angular 
distribution is somewhat shifted toward larger angles as compared to that for charged pions. 
Finally, the distributions on W and arc plotted in Fig. 11 for two E„ intervals. Due to 
a limited resolution of W, as well as the effects of the secondary intranuclear interactions, 
its distribution extends to the region below the nucleon mass. Besides, no peak is seen at 
mjv, due to a small contribution of the identified events of the quasielastic reaction (1) 
as compared to the non-identified one. 
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4 Summary 



The charged current neutrinonuclear reactions ai E„ = 1 - 3 are investigated using the data 
obtained with SKAT propane-frcon bubble chamber irradiated in the neutrino beam at Ser- 
pukhov accelerator. The reconstructed neutrino spectrum is found to be consistent with the 
expected one on the bubble chamber. 

For the first time, in this Ei^ range detailed experimental data are obtained on the character- 
istics of the produced particles. The E^, - dependence of the mean multiplicities of protons, 
charged pions, 7- quanta, 7r° mesons and neutral strange particles and their multiplicity, mo- 
mentum and angular distributions are measured. An indication is obtained, that the yield of 
pions in the neutrinonuclear interactions is reduced as compared to that in i/N interactions, 
which can be qualitatively explaned by their absorption in the nuclear medium. 
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Figure 1: The distribution on: E^^^ for selected events (a), the reconstructed neutrino energy 
(b), the reconstructed energy of the beam neutrino (c). The curve is taken from 
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Figure 2: The distribution for events-candidates for exclusive channels on: the azimutal 
angle between the momenta of the muon and the hadronic system (a); squared transverse 
momentum p'^ (b); ratio (E'"'^ —pl)/EI^^ for events with (j) > 155° and < 0.15 (GeV/c)^ 
(c); El, for exclusive channels (l)-(3) (d); E^, for channels: (1) (full circles), (2) (empty 
circles) and (3) (triangles) (e); invariant hadronic mass W for exclusive channels (l)-(3) 
(full circles) and the channel (2) (empty circles) (f). The curves in Fig.2e correspond to the 
expected dependences for neutrino interactions with (quasi) free nucleons (see the text). 
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Figure 4: The E^, - dependence of the mean momentum of: muons (a); protons (b); charged 
pions (c); 7- quanta (d). 
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Figure 5: The E^, - dependence of the mean energy fracton of the neutrino carried by 
muons and charged pions (a); protons (b); 7- quanta (c); undetected neutral particles (d). 
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Figure 6: The E^, - dependence of the mean values of W (a) and Q"^ (b) 
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Figure 8: The momentum and angular distributions of muons for two Ey intervals. 
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Figure 9: The momentum and angular distributions of charged pions for two E^, intervals. 
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Figure 10: The momentum and angular distributions of 7- quanta corrected (close circles) 
and not corrected (open circles) for their detection efficiency. 
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Figure 11: The distributions on W (a) and (b) 
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